Electronic theory for the normal state spin dynamics in S^RuO^ 
anisotropy due to spin-orbit coupling 
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Using a three-band Hubbard Hamiltonian we calculate within the random-phase-approximation 
the spin susceptibility, \(c\,uj), and NMR spin-lattice relaxation rate, 1/Ti, in the normal state 
of the triplet superconductor Sr2Ru04 and obtain quantitative agreement with experimental data. 
Most importantly, we find that due to spin-orbit coupling the out-of-plane component of the spin 
susceptibility \ zz becomes at low temperatures two times larger than the in-plane one. As a conse- 
quence strong incommensurate antiferromagnetic fluctuations of the quasi-one-dimensional xz- and 
yz-bands point into the z-direction. Our results provide further evidence for the importance of spin 
fluctuations for triplet superconductivity in S1-2RUO4. 
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The spin-triplet superconductivity with X" C =1.5K ob- 
served in layered Sr2Ru04 seems to be a new example 
of unconventional superconductivity Q . The non s-wave 
symmetry of the order parameter is observed in several 
experiments (see for example (2]]j|). Although the struc- 
ture of Sr2RuC>4 is the same as for the high-T c super- 
conductor La2- a; Sr a ;Cu04, its superconducting proper- 
ties resemble those of superfiuid 3 He. Most recently it 
was found that the superconducting order parameter is 
of p- wave type, but contains line nodes half-way between 
the RuC>2-planes [HJ^]. These results support Cooper- 
pairing via spin fluctuations as one of the most proba- 
ble mechanism to explain the triplet superconductivity 
in Sr2RuC>4. Therefore, theoretical and experimental in- 
vestigations of the spin dynamics behavior in the normal 
and superconducting state of Sr 2 Ru04 are needed. 

Recent studies by means of inelastic neutron scatter- 
ing(INS) H and nuclear magnetic resonance (NMR) |7]] 
of the spin dynamics in Sr2Ru04 reveal the presence of 
strong incommensurate fluctuations in the Ru02-plane at 
the antiferromagnetic wave vector Qj = (2tt/3, 2tt/3). As 
it was found in band structure calculations [g] , they result 
from the nesting properties of the quasi-one-dimensional 
d xz - and d yz -bands. The two-dimensional d^-band con- 
tains only weak ferromagnetic fluctuations. The ob- 
servation of the line nodes between the Ru02-plancs 
suggests strong spin fluctuations between the RuC>2- 
planes in z-direction PdTl|. However, inelastic neutron 
scattering observes that magnetic fluctuations are 
purely two-dimensional and originate from the RuC>2- 
plane. Both behaviors could result as a consequence of 
the magnetic anisotropy within the RuCvplane as indeed 
was observed in recent NMR experiments by fshida et 
al. [ p^[ . In particular, analyzing the temperature depen- 
dence of the nuclear spin-lattice relaxation rate on 17 O in 
the Ru02-plane at low temperatures, they have demon- 
strated that the out-of-plane component of the spin sus- 
ceptibility can become almost three time larger than the 



in-plane one. This strong and unexpected anisotropy dis- 
appears with increasing temperature |jl3f . 

In this Communication we analyze the normal state 
spin dynamics of the Sr2Ru04 using the two-dimensional 
three-band Hubbard Hamiltonian for the three bands 
crossing the Fermi level. We calculate the dynami- 
cal spin susceptibility x(q, to) within the random-phase- 
approximation the and show that the observed magnetic 
anisotropy in the Ru02-plane arises mainly due to the 
spin-orbit coupling. Its further enhancement with low- 
ering temperatures is due to the vicinity to a magnetic 
instability. Thus, we demonstrate that as in the super- 
conducting state jl4| the spin-orbit coupling plays an im- 
portant role also for the normal state spin dynamics of 
Sr 2 Ru04. We also discuss briefly the consequences of 
this magnetic anisotropy for Cooper-pairing due to the 
exchange of spin fluctuations. 

We start from the two-dimensional three-band Hub- 
bard Hamiltonian 

H = H t + H v = J~] frda^; .ak,;<7 + ^ Uin^rim, 

k,<7 / i,l 

(1) 

where ak,z CT is the Fourier-transformed annihilation oper- 
ator for the di orbital electrons (I = xy, yz, zx) and Ui 
is the corresponding on-site Coulomb repulsion. t\a de- 
notes the energy dispersions of the tight-bindings bands 
calculated as follows: tyi = — eo — 2t x cos k x — 2t y cos k y + 
it' cos k x cos k y . We choose the values for the parame- 
ter set (e ,t x ,t y ,t') as (0.5, 0.42, 0.44, 0.14), (0.24, 0.31, 
0.045, 0.01), and (0.24, 0.045, 0.35, 0.01)eV for d xy -, d zx -, 
and dj,2-orbitals in accordance with band-structure cal- 
culations |l5| . The electronic properties of this model in 
application to Sr2Ru04 were studied recently and as was 
found can explain some features of the spin excitation 
spectrum in Sr2Ru04 P Jl4| , |T6| , |lT| . However, this model 
fails to explain the observed magnetic anisotropy at low 
temperatures Jli|] and line nodes in the superconducting 
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FIG. 1. Calculated Fermi surface for a RuC>2 plane in 
Sr2Ru04 taking into account spin-orbit coupling. 

order parameter below T c which are between the Ru02- 
planes. On the other hand, it is known that the spin-orbit 
coupling plays an important role in the superconducting 
state of in S^RuCU Jlij . This is further confirmed by 
the recent observation of the large spin-orbit coupling in 
the insulating Ca2RuC>4 0. Therefore, we include in 
our model spin-orbit coupling: 



A L^Si 



(2) 



where the angular momentum Lj operates on the three 
^q-orbitals on the site i. Similar to an earlier approach 
{[jj l, we restrict ourselves to the three orbitals, ignoring 
e2 9 -orbitals and choose the coupling constant A such that 
the t2 9 -states behave like an I = 1 angular momentum 
representation. Moreover, it is known that the quasi- 
two-dimensional xy-band is separated from the quasi- 
one-dimensional xz- and j/z-bands. Then, one expects 
that the effect of spin-orbit coupling is small and can be 
excluded for simplicity. Therefore, we consider the ef- 
fect of the spin-orbit coupling on xz- and yz-bands only. 
Then, the kinetic part of the Hamiltonian H t + H so can 
be diagonalized and the new energy dispersions are 



e k,yz = (*k,j/^ + tk,xz + Ak)/2 
€ k,xz = (^k,t/2 + tk.xz ~ ^-k)/2 



(3) 



where — yj (t^ yz — t^ xz ) 2 + A 2 , and a refers to spin 
projection. One clearly sees that the spin-orbit coupling 
does not remove the Kramers degeneracy of the spins. 
Therefore, the resultant Fermi surface is consists of three 
sheets like observed in the experiment. Most impor- 
tantly, spin-orbit coupling together with Eq. (1) leads to 
a new quasiparticle which we label by pseudo-spin and 
pseudo-orbital indices. The unitary transformation 
connecting old and new quasiparticles is defined for each 
wave vector and lead to the following relation between 
them: 
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FIG. 2. Diagrammatic representation for the transverse 
and longitudinal components of the magnetic susceptibility. 
The full lines represent the electron Green's function with 
pseudospin a and pseudo-orbital Z-indexes. g+ and g z denote 
the vertexes as described in the text. 
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where u m k 



and VmU 



„T4 



\/(tk,yz-t k , X zTAk) 2 +\ 2 

The '-' and '+' signs refer to the 



772 = 1 and to = 2, respectively. 

In Fig.l we show the resultant Fermi surfaces for each 
obtained band where we have chosen A = lOOmeV in 
agreement with earlier estimations p^ , p^7| . One immedi- 
ately sees that xz- and yz-bands split around the nested 
parts in good agreement with experiment pq |. Thus, 
spin-orbit coupling acts as a hybridization between these 
bands. However, in contrast to hybridization spin-orbit 
coupling introduces also an anisotropy for the states with 
pseudo-spins t and J,. This will be reflected in the mag- 
netic susceptibility. Since the spin and orbital degrees 
of freedom are now mixed in some spin-orbital liquid, 
the magnetic susceptibility involves also the orbital mag- 
netism which is very anisotropic. 

For the calculation of the transverse, xl > an d longitu- 
dinal, x\ z i components of the spin susceptibility of each 
band I we use the diagrammatic representation shown in 
Fig. 2. Since the Kramers degeneracy is not removed by 
the spin-orbit coupling, the main anisotropy arises from 
the calculations of the anisotropic vertex g z = l z + 2,s 2 
and g + = l + + 2s + calculated on the basis of the new 
quasiparticle states. In addition, due to the hybridization 
between xz- and yz-bands we also calculate the trans- 
verse and longitudinal components of the the interband 
susceptibility xw ■ Then, for example, 
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FIG. 3. Results for the real part of the out-of-plane (solid 
curve) and in-plane (dashed curve) magnetic susceptibilities, 
Re x(q, w), calculated within RPA using U = 0.575eV along 
the route (0, 0) -> (tt, 0) -> (tt, it) -> (0, 0) within the first 
Brillouin Zone at temperature T = 100K. 



M2kM2k+q + V2k«2k+q + V2(u2k«2k+q + «2kU2k+q) 



£ kx Z - £ k+q^ + W + l0 + 



(6) 



where /(cc) is the Fermi function and ti£ and v%. are 
the corresponding coherence factors which we have 
calculating through the corresponding vertexes using Eq. 
(4). For all other orbitals the calculations are straight- 
forward. Note, that the magnetic response of the xy- 
band remains isotropic. One clearly sees the difference 
between longitudinal and transverse components which 
results from the calculated matrix elements. Moreover, 
the longitudinal gets an extra term due to l z while the 
transverse does not contain the contributions from l + or 
The latter occur due to the fact that xz- and yz-states 
are a combination of the real orbital states |2, +1 > and 
1 2,-1 >. Thus the transition between these two states 
are not possible with l + or i_ operators. Therefore, each 
component of the longitudinal susceptibility gets an ex- 
tra term in the matrix element that sufficiently enhances 
their absolute values. 

Assuming C/y = SijU one gets the following expres- 
sions for the transverse susceptibility within RPA: 



xti (q> w ) 



1 - Uxo t i (<l' u ) 
and for the longitudinal susceptibility 

XrpaM^) = 

1 -u 2 xi.M^)xlMiu) 



(7) 



(8) 



In Fig. 3 we show the results for the real part of the 
transverse and longitudinal total susceptibility, Xrp'a ~ 
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FIG. 4. Temperature dependence of the imaginary part 
of the spin susceptibility divided by uj s f and summed over q. 
Note, zz and +— refer to the out-of-plane (solid curve) and 
in-plane (dashed curve) components of the RPA spin suscep- 
tibility. In the inset we show the corresponding frequency 
dependence of the ImxflPA(Qi, w) at the IAF wave vector 
Qi = (27r/3, 2n/3). The results for the out-plane component 
(solid curve) are in a quantitative agreement with INS exper- 
iments 



Y^iXRPAi alon g the route (°>°) -> K°) {k,-k) -> 

(0, 0) in the first Brillouin Zone for U = 0.505eV . Note, 
the important difference between the two components. 
The longitudinal component of the spin susceptibility is 
almost three times larger than the transverse one all over 
the Brillouin Zone. Moreover, despite of some structure 
seen in Xrpa a ^ (27i"/3, 27r/3) there are no real incommen- 
surate antifcrromagnetic fluctuations at this wave vec- 
tor. On the other hand, the structure in Xrpa a ^ the 
same wave vector refers to real fluctuations. The lat- 
ter is seen in the inset of Fig. 4 where we present the 
results for the frequency dependence of the imaginary 
part of the total susceptibilities at Q; = (2-7r/3, 27r/3) 
and temperature T=20K. The longitudinal component 
reveals a peak at approximately u s f = 6meV in quanti- 
tative agreement with experimental data on INS |(| . The 
transverse component is featureless showing the absence 
of the incommensurate antifcrromagnetic spin fluctua- 
tions. Thus, the fluctuations in the transverse suscep- 
tibility are isotropic and ferromagnetic-like. Therefore, 
antiferromagnetic fluctuations are present only perpen- 
dicular to the Ru02-plane. 

We also note that our results are in accordance with 
earlier estimations made by Ng and Sigrist [jl9| with one 
important difference. In their work it was found that 
the IAF are slightly enhanced in the longitudinal compo- 
nents of the xz- and yz-bands in comparison to the trans- 
verse one. In our case we have found that the longitu- 
dinal component of the magnetic susceptibility strongly 
enhances due to otbital contributions. Moreover, we 
show by taking into account the correlation effects within 
random-phase-approximation(RPA) the IAF are further 
enhanced in the z-direction. 
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FIG. 5. Calculated normal state temperature depen- 
dence of the nuclear spin-lattice relaxation rate T^ 1 of 17 in 
the Ru02-plane for the external magnetic field applied along 
c-axis (dashed curve) and along the ab- plane (solid curve). 
Down- and up-triangles are experimental points taken from 
Ref. |l3| for the corresponding magnetic field direction. 



In order to see the temperature dependence of the mag- 
netic anisotropy induced by the spin-orbit coupling we 
display in Fig. 4 the temperature dependence of the quan- 



tity E q 1 



™-XRPA(<j,U s f) 



for both components. At room 
temperatures both longitudinal and transverse suscepti- 
bilities are almost identical, since thermal effects wash 
out the influence of the spin-orbit interaction. With de- 
creasing temperature the magnetic anisotropy arises and 
at low temperatures we find the important result that 
the out-of-plane component x zz is about two times larger 
than the in-plane one (x zz > X H /2). 

Finally, in order to compare our results with experi- 
mental data we calculate the nuclear spin-lattice relax- 
ation rate for 1 O ion in the RuC>2-plane for different 
external magnetic field orientation (i = a, 6, and c) 



1 



W 2 - 1 ql 



(9) 



where is the g-dependent hyperfine-coupling constant 
perpendicular to the z-direction. 

In Fig. 5 we show the calculated temperature depen- 
dence of the spin-lattice relaxation for an external mag- 
netic field within and perpendicular to the Ru02-plane 
together with experimental data. At T — 250K the spin- 
lattice relaxation rate is almost isotropic. Due to the 
anisotropy in the spin susceptibilities arising from spin- 
orbit coupling the relaxation rates become different with 
decreasing temperature. The largest anisotropy occurs 
close to the superconducting transition temperature in 
good agreement with experimental data p3[ . 

To summarize, our results clearly demonstrate the es- 
sential significance of spin-orbit coupling for the spin- 
dynamics already in the normal state of the triplet super- 
conductor Sr2RuC>4. We find that the magnetic response 
becomes strongly anisotropic even within a Ru02-plane: 



while the in-plane response is mainly ferromagnetic, the 
out-of-plane response is antiferromagnetic-like. 

Let us also remark on the implication of our results 
for the triplet superconductivity in Sr 2 Ru04. In a previ- 
ous study |ll]] , neglecting spin-orbit coupling but includ- 
ing the hybridization between xy- and xz—, j/z-bands, we 
have found ferromagnetic and IAF fluctuations within 
the a6-plane.This would lead to nodes within the RuC>2- 
plane. However, due to the magnetic anisotropy in- 
duced by spin-orbit coupling, a nodeless p— wave pair- 
ing is possible in the RuC>2-plane as experimentally ob- 
served. Our results provide further evidence for the im- 
portance of spin fluctuations for triplet superconductivity 
in Sr 2 RuC>4. 
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